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I/"") ' Abstract 

We investigate the behavior of entanglement-entropy on a broad 
f--*) ' scale, that is, a large class of systems, Hamiltonians and states de- 

. scribing the interaction of many degrees of freedom. It is one of our 

T 7*^ ' aims to show which general characteristics are responsible for the dif- 

JL . ferent types of quantitative behavior of entantglement-entropy. Our 

main lesson is that what really matters is the degree of degeneracy 
C"| , of the spectrum of certain nearby reference Hamiltonians. For calcu- 

lational convenience we study primarily systems defined on large but 
• rH , finite regions of regular lattices. We show that general vector states, 

f\ ' being not related to some short-range Hamiltonian do not lead in the 

generic case to an area-like behavior of entanglement-entropy. The sit- 
uation changes if eigenstates of a Hamiltonian with short-range interac- 
tions are studied. We found three broad classes of eigenstates. Global 
groundstates typically lead to entanglement-entropies of subvolumcs 
proportional to the area of the dividing surface. Macroscopically ex- 
cited (vector)states have in the generic case an entanglement-entropy 
which is proportional to the enclosed subvolume and, furthermore, dis- 
play a certain Gibbsian behavior. Low-lying excited states, on the 
other hand, lead to an entanglement-entropy which is usually propor- 
tional to the logarithm of the enclosed subvolume times the area of 
the dividing surface . Our analysis is mainly based on a combination 
of concepts taken from the perturbation theory of Hamiltonians and 
certain insights coming from the foundations of quantum statistical 
mechanics. 



1 Introduction 



The microscopic origin of entropy in, for example, black hole physics and 
the so-called area law is still to some extent kind of a mystery. Or more 
precisely, there exist a variety of different explanations. For a nice discussion 
in form of a Galilean trialogue see pp. Another very readable overview is 
given in [21], a collection of theses can also be found in [3]. Of perhaps even 
more interest is the proposed holographic principle (and the range of its 
validity), which discusses the emergence of entropy on a very fundamental 
level of physics, that is, "empty" space (-time) and its vacuum fluctuations. 
We do not attempt to give a complete list of references (as we focus in the 
following on a related but slightly different question), cf. for example the 
nice review [3] and the many references therein. 

What is at least clear is that the black hole horizon divides space-time 
into exterior and interior regions. Hence, one might venture the idea that 
a basic role in this question should be played by some form of entanglement 
of two (in some macroscopic respect) separated regions. Such a form of 
entanglement can however be realized on different levels of our theoretical 
description, ranging from more ordinary ones (standard quantum mechanics 
or quantum field theory) to more pristine ones ( e.g. the notorious Planck 
scale). On this point people widely disagree at the moment. As to our 
personal point of view, we think that the holographic principle indicates a 
drastic change of the statistical mechanical preassumptions usually taken 
for granted on the more ordinary scales as e.g. locality, good clustering con- 
ditions or sufficient decay of correlations etc., if one enters the fundamental 
regime of the Planck scale. Some very sketchy remarks can be found at the 
end of [5], a more detailed analysis is forthcoming (|6j). Be that as it may, 
one should emphasize the following: 

Observation 1.1 The following two situations are different from a logical 
point of view, for one, entanglement- entropy induced by a fixed global state 
as e.g. the groundstate of a given Hamiltonian, for another, the maximally 
possible quantum of entropy or information, S max , which can be stored in 
a given volume, the question which is adressed in the holographic principle. 
In the latter case also the highly excited states of a Hamiltonian should be 
investigated, a problem we explicitly address in the following sections. 

In any case, in a first step, it seems to be useful to investigate ordinary 
models of quantum (field) theory and see to what extent their entropic 
behavior reflects some of the properties condensed in notions like area law 
and/or holographic principle. A certain kind of area-law- like behavior shows 
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up already in ordinary quantum field theory, see e.g. [7] or [8]. Similar 
phenomena were also observed in [9], [TO] and more recently in [5]. As this 
is an interesting point, we will treat it in more detail elsewhere. 

To make the problem more accessible, we start with the investigation of 
systems of many degrees of freedom living on a large regular lattice (of ar- 
bitrarily small but finite lattice constant). In [11] or [12] or the more recent 
[13j large arrays of coupled harmonic oscillators and their respective ground- 
states have been studied, with some sort of continuum limit performed in 
the end, yielding the vacuum state of e.g. a Klein-Gordon field theory. Note 
however that taking a continuum limit is quite delicate since, in order to 
get a finite result, some form of cut-off is necessary (as to other possible 
interesting applications to quantum field theory see for example [14J). 

As everybody knows, the harmonic oscillator is an extremely well-behaved 
quantum system. The same holds for regular arrays of coupled harmonic 
oscillators. There exists a large arsenal of methods to extract useful in- 
formation from these model systems. But nevertheless, the question of 
entanglement- entropy of e.g. the groundstate of this well-behaved model 
system turns already out to be quite intricate. This may, among other 
things, have its roots in the observation that the concept of entanglement- 
entropy and its quantitative behavior is not so easy to visualize, one reason 
being, that it is not a truely local concept. 

It was found that if one divides a large volume, V , into V\ U Vi and 
restricts the groundstate, over V to the subvolumes Vi,V2, the cor- 
responding density matrices, W\,W2, have an entropy, S = — ^Wjlntfj, 
which is proportional to the area of the dividing surface between V\ and 
V2I One may therefore speculate that this observation has something to do 
with the area law in black hole physics and corresponding arguments were 
advanced in the two above cited papers. 

If one persues such an idea, various questions immediately suggest them- 
selves. 

i) In what respect is the groundstate of the Hamiltonian of a large system 
(i.e., many degrees of freedom being involved) a particular state in this spe- 
cial context. Put differently, what particular property of a global vectorstate 
(pure state) influences the entropic behavior of the partial traces, W\, W2? 

ii) In what respect does the picture change if we go over to lowly or highly 
excited eigenstates of such a Hamiltonian? 

iii) Is the result, found for an array of harmonic oscillators, generic, that 
is, is it to a high degree independent of the particular model system being 
studied and what are the really important prerequisites? 

These are the questions we will adress in the following sections. A prelimi- 
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nary qualitative answer may perhaps already be given. Perhaps a little bit 
surprisingly, what really seems to matter is the degree of degeneracy of the 
spectral values of certain nearby reference Hamiltonians. For groundstates, 
for example, this means that the area-behavior is a result of their usual 
non- degeneracy. 

We begin our investigation with defining the general framework of sys- 
tems, living on large regular lattices. We then discuss in section 3 several 
general properties of the notion of entanglement- entropy. We mention for 
example that a Gibbs state over V\ can always be extended to a pure vector 
state, ^, on V by enlarging the system on V% with V = V\ U V 2 . That is, the 
Gibbs state is now the partial trace of a global vectorstate. As, typically, 
such Gibbs states have an entropy which is proportional to the volume, V\, 
this shows that, in general, entanglement-entropy can not be proportional 
to the area of the surface, dividing V\ and V 2 - We learn from this simple 
observation that, in order to get such a special result, we do need some more 
specifications. 

This observation teaches us yet another lesson. One sometimes hears the 
argument that entanglement-entropy has to be proportional to the area of 
the dividing surface as this is the only quantity which is the same for both 
subvolumes, V\ and V 2 while the volumes themselves can of course differ 
from each other (note that, perhaps a little bit surprisingly, entanglement- 
entropy is always the same for Vi and V 2 ). This argument is obviously 
incorrect. What in general may happen is that we have something like the 
following 

S{W 1 ) = a l -\V l \=a 2 -\V 2 \= S{W 2 ) (1) 

with W±,W 2 the density matrices representing the total state, ^, on V\,V 2 
respectively. That is, both entropies may happen to depend linearly on the 
volume but with certain (volume-dependent) adjusting prefactors. As long 
as the volumes are finite, there is no unique answer as to the geometric 
dependencies. We need some kind of thermodynamic limit. In [13] one kind 
of such limit is discussed. One can, for example, keep the dividing surface 
fixed and let the subvolumes approach infinity in the orthogonal direction, 
or, on the other hand, it may become necessary, to scale them in a prescribed 
way with the prefactors remaining basically constant. In the former case we 
encounter a linear dependence on the area of the dividing surface if the 
limit is finite, in the latter (scaling) case we may get a linear asymptotic 
dependence of S(Wi) on the volume. 

At the end of that section we introduce a certain coarse-grained version 
or approximation of the notion of entanglement-entropy for the subvolumes 
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V\,V2, which, in some sense, resembles kind of a microcanonical ensemble 
version of entropy. In our particular case it is related to the logarithm of 
the dimension of certain Hilbert subspaces onto which the pure vectorstate 
happens to be essentially projected when calculating the relative trace. This 
is an important simplification: While the original standard notion of entropy, 
S := — w i ' l n w i > is deceptively simple at first glance, this is a very 
missleading impression. Even in the few cases where the Boltzmann weights, 
Wi, are known, what is usually more important is their (frequently huge) 
degeneracy and/or the local density of states, which in general are not at 
our disposal. 

Remark:We recently learned that this approximation scheme we will discuss 
below has a certain resemblance to an approximation method used in the 
density matrix renormalisation group, DMRG and which can roughly be 
described as finding the essential subspace. It is also of relevance in data 
compressing and the like (see e.g. [31 j ,[32j,[33j,[34j). 

While we attempt to treat the problem of entanglement-entropy on quite 
a broad scale (as far as space dimension and class of hamiltonians are con- 
cerned), there are interesting topics we do not! address. In our scenario we 
deal with a spatial tensor product structure of the underlying Hilbert space. 
Most notably in condensed matter physics there exist applications where 
other kinds of configurations prevail. One example is the case of fermionic 
quantum particles where the antisymmetrization postulate prevents such a 
representation. In that situation one can for example work in a mode repre- 
sentation; (basically Fock space) and discuss entanglement in such a mode 
space. Some work in this direction may be found in the review article [15] . 
[16], [T7] or [18]. We think that parts of our strategy can also be applied in 
this field of research but this is not done in the present paper. 

Recently, Shi ([H]) attempted to extend such a mode- version of entan- 
glement to relativistic quantum field theory. This is certainly an interesting 
programme but we would like to hint at a number of impending difficulties 
and obstacles. In free field theories and Fourier space, occupation num- 
ber representations for example have a well-defined meaning. In coordinate 
space and, a fortiori, interacting theories such an approach is met with great 
difficulties due to the intricacies of the relativistic quantum regime. For one, 
expressions like 4>\x)(j){x) are no longer operators but operator- valued dis- 
tributions. The same applies to eigenstates, "localized at a;". What one can 
use instead are localized algebras of observables living in regions O of space- 
time. One may then restrict e.g. the vacuum state to such local algebras, 
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thus getting a mixed state. 

It is however obvious that in general a naive splitting into, say, a bipar- 
tite system, having kind of a tensor-product structure, does not work (due 
to quantum non-locality and other problems; a catchword being the Reeh- 
Schlieder Property). In any case, the mathematics becomes extremely in- 
volved as most of the naive expressions become infinite (for example traces). 
Only for a few scenarios one can say more (like e.g. the Rindler-wedge) . The 
interested reader may consult section V,5 of the book by Haag ([2Q]),[l3] or 
|21j for a general orientation. 

Furthermore, our arguments rely to a certain extent (at least in the 
form presented here) on short-range correlation and short-range interac- 
tions. This is the field where one may expect that results are generic 
to some extent. Therefore our arguments do not (immediately) apply to 
quantum critical systems, i.e. systems displaying quantum phase transitions 
at temperature zero. These are presently mostly studied in low dimen- 
sions (space dimension frequently being one) by employing usually quite 
model-specific techniques. An incomplete selection of papers is [22] . |23j . [21] . 
[25], [26], [27], [28]. 

Given the generality of our admitted class of models, it is not our aim 
and out of the question to really compute in every case the complete set of 
eigenvalues and their degeneracy or spectral density. This holds the more 
so for the v. Neumann entropy. We rather view the merit of our approach 
in providing a general scheme which only relies on a few characteristics 
of the respective models. As a consequence, our arguments and estimates 
do not yield precise numerical results but typically make statements about 
entanglement-entropy in some leading order (for more details see below). 

For the same reason (a large class of different admissible Hamiltonians) 
we develop in section d] in quite some detail various concepts belonging to 
the perturbational analysis of Hamiltonians. This turns out to be a very 
intricate subject matter in our context. The main reason is that, due to the 
dense distribution of eigenvalues typical for Hamiltonians of a large number 
of degrees of freedom, ordinary perturbation theory is only of a limited value. 

In the last three sections we then deal with our main topic in a more 
quantitative way. In section [5] we treat the case of the groundstate of a 
Hamiltonian with finite range interactions. In the following section we anal- 
yse the situation for eigenstates which are highly excited, that is, for energies 
which are a macroscopic distance apart from the groundstate energy. In the 
last section we study, on the other hand, eigenvalues of the Hamiltonian 
which lie in the vicinity of the groundstate energy. We find the following: 
i) Groundstates come with an entanglement-entropy which is proportional 
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to the area of the dividing surface. 

ii) Highly (i.e. macroscopically) excited states lead in the generic case to 
entanglement-entropies which are proportional to the volume of the system. 
The notion generic is explained in the respective section. Furthermore, the 
restriction of the pure vectorstate to the respective subvolumes share a cou- 
ple of properties with ordinary temperature states. As for more detailed 
qualifications see the corresponding section. 

iii) Low-lying excited eigenstates states have an entanglement-entropy which 
is roughly proportional to the logarithm of the corresponding volume times 
the area of the dividing surface. 

To sum up what we think have been important ingredients in our analy- 
sis, the arguments strongly rely on the finite range of the interactions, that 
is some version of locality. Furthermore, we assume the system to be away 
from a (quantum-) critical state where local degrees of freedom are typically 
long-range correlated. In this special regime results happen to be somewhat 
model dependent and our more general arguments need not automatically 
apply. As to the applied methods, two tools stand out. First, a combi- 
nation of various perturbational arguments, second, a couple of technical 
ideas which have proved their worth already in the foundations of statistical 
mechanics. 

We want to add a further remark. One might perhaps get the idea that 
high entanglement is necessarily related to long-range correlations. That 
this is not the case can be inferred from our following investigations. It is 
for example easy to construct a highly entangled eigenstate of a short-range 
Hamiltonian. What is only important is that the respective eigenvalue hap- 
pens to be highly degenerate or that the full Hamiltonian is the perturbation 
of such a Hamiltonian with highly degenerate eigenstates. This is exploited 
in section [6j Note that, typically, most eigenvalues of Hamiltonians of many 
degrees of freedom are highly degenerate, an exception being usually the 
ground state. For general vector states there is even less correlation be- 
tween entanglement and range of correlations. An interesting observation 
in this direction is for example also made in [29J and [30] where these two 
notions are compared with each other. 

A related line of ideas is pursued in [35] . Again Shi concentrates mainly 
on mode expansions of many body systems and studies in particular the 
regime of phase transitions and spontaneous symmetry breaking. For exam- 
ple, in the case of Bose-Einstein condensation he remarks that a macroscopic 
occupation of the ground state leads to what he calls disentanglement with 
the rest of the modes. This is reasonable as the ground state occupation 
number operator becomes almost a classical observable in this situation. He 



6 



further addresses the relation between entanglement and range of correla- 
tions in this context. It conforms with our observations that long-range 
correlations make it in general difficult to divide a system into, for exam- 
ple, a weakly correlated bipartite system, while high entanglement does not 
necessarily imply long-range correlations. 

2 Some General Framework and Estimates 

Concerning this section a standard reference is for example [36]. To avoid 
irrelevant complications it is reasonable to choose as the general context dis- 
crete dynamical models living on a regular lattice, ZJ^, with lattice constant 
a and space dimension d. The points of the lattice are denoted by Xj, the 
Hilbert space sitting at Xj by Tii. Tii may be finite or infinite dimensional 
(as is for example the case for an array of harmonic oscillators). If all the 
Tii have the same finite dimension it is denoted by D. 

In the following we are interested in large but finite subsets of Z^, denoted 
by V. For convenience V may be chosen as a (i-dimensional cube of side 
length L, containing Ny = (L/a) d lattice points. The general case of regions 
of arbitrary shape can of course be treated in more or less the same way 
apart from some (irrelevant) numerical and technical details. The Hilbert 
space, TLy, over V is the tensor product 

■Hv=®rLi (2) 

spanned by tensor-monomials 

v = vi <g> • • • <g> vn v , Vi eHi (3) 

If one chooses a basis, {ej*}, in each Hi, general vectors in TLy are linear 
combinations of the basic vectors 

4 1 <g> • • • ® e^; , 1 < k < D (4) 

(in case the dimension of all Tii is finite and equal) . Then the dimension of 
TLy is D Nv . The scalar product on TLy is given by 

(«l <g> • • • ® vavK ® • • • ® v' Nv ) := IJ(uiK) (5) 
and (multi) linearly extended. 
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Hamiltonians are now denned as follows. If A4 are operators (or matrices) 
on Hi, operator monomials on TLy are given by 

(gMi 0^o Vi (6) 

Xi 

and (multi) linearly extended on more general vectors. Refraining from a 
too abstract approach we restrict ourselves to Hamiltonians given by a sum 
over one-point, two-point,..., k-point-interactions. This means 

Xyo*). j> 2 (*,, *,),••• (7) 

where 

s 1 ^) = - * 2 (^) = A (2) ® Af , ... (8) 

Assumption 2.1 To avoid unneccessary complications we assume that only 
interaction terms up to a certain finite order occur in the Hamiltonian and 
that all the occurring operators A{ are hermitean. 

As all the Ai,Aj do commute for i 7^ j, the interaction terms defined 
above are hermitean as is their sum. Typical examples are spin systems. 
One of the simplest examples is 

H = J ■ ^2 ~Oi ® o^i+i) =: J • ~a i ~a { i +1 ) (9) 

i i 

with a denoting the usual Pauli-matrices and extension to higher dimensions 
being straightforward. 

Usually one makes the additional assumption that the interaction has a 
finite range and is translation invariant. 

Definition 2.2 The interaction is called translation invariant if with 
<& fc (xi 1 , . . . , Xi k ) also $ fc (xj x + r, . . . , Xi k + r) occurs in the Hamiltonian with 
r some lattice vector and both being the same operators. 

Definition 2.3 The interaction is called to have a finite range, 5, if with 
one lattice point fixed, e.g. x,^, only finitely many members ^(x^, yj 2 , ... ,yi k ) 
different from zero exist in H when the y^ vary over Zf , all yj having a dis- 
tance from Xi which is < 5. 
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We not that we give these details only for consistency reasons, whereas 
only the general properties of such Hamiltonians will be of relevance in the 
following like their translation invariance and the finite range. 

What will be of interest for the following discussion is the restriction of a 
general Hamiltonian, H, to a certain subvolume V C Z„ or the correspond- 
ing sub-Hilbert-space Hy. 

Definition 2.4 The restriction of a general Hamiltonian, H, of the kind 
described above to a subvolume V C 2^ or Tiy is the operator Hy which con- 
sists of all the interaction terms which lie in V , i.e., only k-tuples (x^ , . . . , Xj fc ) 
are admitted with all the x^ £ V. 

Hy.= Y, *W..,s ifc )=:$> fc ( 10 ) 

k;x it £V n(V) 

(n(V), the number of terms in the sum being of order 0(\V\) ). 

Another important part of the total H relative to a given volume V is the 
boundary contribution Hy 1 . 

Definition 2.5 By Hy - we denote the part of the total H which consists of 
all interaction terms which have both lattice points in V and the dual set 
V':=Z d a \V. 

H b v d := £ ^ k {x ll ,...;y Hl ...,y %k ) (11) 

k;xeV,yeV 

Observation 2.6 It will be important in the following that with the Hamil- 
tonian having finite range and k < K, the number of terms in H'y is of 
order (area of boundary ofV). We denote this number of terms by n(V, V). 

We can make more detailed statements if we concentrate on the large sub- 
class of models with all the $ fc bounded operators (which holds for example 
if all the Tii are finite dimensional). This together with the assumed trans- 
lation invariance and the above observation allow us to make the following 
important operator estimates. 

Proposition 2.7 Under the preceding assumptions we have 

|| Hy || < C ■ n{V) , || Hy \\<C ■ n(V, V') (12) 

where in the generic case, which we usually assume to prevail, both norms 
are actually proportional to the volume, the area of the boundary, respec- 
tively. 
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In case the local Hilbert spaces Hi are finite dimensional, all the occurring 
Hamiltonians are large but nevertheless finite hermitean matrices for finite 
V , hence having a discrete spectrum, with the number of eigenvalues (count- 
ing multiplicity) being the same as the dimension of the respective Hilbert 
spaces (the eigenvalues being the zeros of the characteristic polynomial). 

It is obvious that this latter observation greatly simplifies the quantitative 
analysis. 

We will in the following mainly discuss the case where the local Hilbert 
spaces have finite dimension. As the extension of the results to the case 
where the local Hilbert spaces, Hi have infinite dimension is not entirely 
straightforward (like e.g. arrays of coupled harmonic oscillators), we make 
some corresponding remarks in appendix [Al 

3 Some Properties of the Partial Trace and En- 
tanglement Entropy 

Let us take a vector state, from the Hilbert space Tty with V large (or 
macroscopic). Let us divide V into V1UV2 with the respective Hilbert spaces, 
H±,H 2 - We have 



which defines the states u\ , uj 2 on 7ii,TC 2 . Choosing a basis e« (g) e'- in the 
tensor product Tiy = TL\ <8> TL 2 we get 



H v = Hi ® U 2 , Ay = Ai <8> A 2 



(13) 



We can evaluate the pure state, $ on the restrictions 7Yi,7^2 : 



:= ($|(^1 (1) ® 1)$) , u 2 {B {2) ) := ($|(1 <g> 5( 2 ))$) (14) 




(15) 



and 



1)*) 




y~]bji(ei\Aej) 



(16) 



with 



dim(H2) 




(17) 



m=l 



This can be rewritten as 




(18) 
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The following result is remarkable and leads to the definition of the 
notion of entanglement- entropy. 

Theorem 3.1 The partial traces W\ , W2 have the same spectrum and the 
eigenvalues 7^ have the same multiplicity while the respective zero-eigenvalues 
may have different degeneracies, depending on the in general different di- 
mensions of Hi , H2- 

Corollary 3.2 The (entanglement-) entropies of the states u>i , UJ2 orWi , W2 
are the same, i.e. 

S 2 (W 2 ) = Si(Wi):=-^\ i -ln\ i (19) 

while the respective volumes or dimension of Hilbert spaces can be very dif- 
ferent. 

As the full proof of this important result is frequently omitted in the liter- 
ature or is incomplete we give, for the convenience of the reader, our own 
version in appendix [Bl 

On the other hand we know that the usual thermodynamic entropy is 
an extensive quantity and depends in general linearly on the volume. So let 
us assume we have a Gibbs-state on the volume V\ with Hamiltonian H\, 
its eigenvalues and eigenstates being Ei , ^i. That means: 

LOi(A) := Tr{e- pHl ■ A)/Tr(e^ Hl ) (20) 

with (3 the inverse temperature 1/kT. It follows 

ft = - 5> , Pi = e-^/E e ~ PEl ( 21 ) 

and 

Si(Vi) ~ Vi (22) 

We now adjoin a volume V2 with a Hilbert space H2 of sufficiently high 
dimension. We define the following vector, ^, in Hi <S> Hi- 

* -=^y/Pi^i®ei (23) 
{ei spanning a basis in H2 and degeneracies being included). We have 

= Tr{e- pHl ■ A^)/Tr{e- pHl ) (24) 

We conclude: 
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Observation 3.3 Every Gibbs state over Ai can be represented by a vector 
state in a sufficiently large ambient Hilbert space. Restricted to Hi this 
vector state is the partial trace, but we see from our above reasoning that the 
entanglement- entropy of ^ relative to V\ is now proportional to the volume 

It is important for the physical understanding and intuition to get a 
better feeling how entanglement-entropy is affected in both a quantitative 
and qualitative way by different physical conditions. It seems to be useful 
to introduce concepts which are, at least in a rough sense, of a similar 
character as entanglement-entropy because the latter can only be calculated 
in very few special cases. We think in this context of the relation between, 
for example, the canonical and the microcanonical ensemble in statistical 
mechanics. 

For one, it is obvious that the entropy of a state, ^o, reduced to, say, V\ or 
V2, does not depend on the choice of Hilbert space bases in the corresponding 
Hilbert spaces. So, in the following, we will frequently subdivide V\ or V2 
further into, say, V\ = V( U V" with V" a boundary layer in V\ neighboring 
upon the interface, separating V\ and V2. One can then equally well choose 
a basis in Hi by forming the tensor product of the respective bases in H'i 
and H'[ etc. 

For another, a famous theorem of E.Schmidt and von Neumann (|41|.|42|) 
states that ^0 in Hi <g> H2 can be represented in the special form 

*o = X}v / V0i,i®0i,2 (25) 

i 

with fc i , (j>i 2 particular orthonormal bases in Hi , H2 (the so-called Schmidt- 
basis; this was also exploited in |43j). 

Remark: In modern parlance this is nothing but the theorem that a com- 
pact operator can be put into such a canonical spectral form. Note in this 
context that ^q, viewed as an operator from H2 to Hi, belongs to the 
Hilbert- Schmidt- class as a consequence of its normalisation as a vector. 

^0 reduced to Hi then yields: 

(*o\AW* ) = *i ■ = Tr(Wi ■ A) (26) 

i 

What is important in this particular representation is that both systems 
of vectors are orthonormal. With $ = Cijei®e'j we can of course always 
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write 



(27) 



with 



52^4 • r > 52 



(28) 




or 



$ = 52 II II ' ei ® & = 52 H ^' 



(29) 




(f>i,ipj being the corresponding unit vectors. 
The reduction to TL\ , TL2 then yields: 



(*|A( 1 )*) = ^(v J -iM0. (*|£ (2) $) 



52(&i B &) 



(30) 



or 



($i^w$) = 5] 11 ^ 11 2 •(^•im-) , ($ib( 2 )$) = x; 11 ^ 11 2 



Remark:Note that the state, given as an expectation functional, is of course 
identical to the state in e.g. formula (|18p . But now the way the state is 
concentrated in a subspace (either spanned by the tpj or the (pi) of the total 
space is made transparent. Note that the subspace may depend on the 
chosen basis in contrast to its dimension (see below) 

Trying to relate the canonical version of the notion of entropy, S(W) := 
— ^2iPi Inpi with some other (perhaps coarser) concept, we were inspired by 
the beautiful analysis of the entropy concept, as it is laid out in [44], sect. 7 of 
chapt.l. It is shown there that the natural quantity which is relevant in this 
context is the number of microscopic quantum states, Ar, a macrostate is 
smeared over or, in other words, the number of microstates which essentially 
contribute in a macrostate (taken with equal weights). It is then shown in 
[44j that in the regime of equilibrium statistical mechanics, the logarithm of 
this quantity coincides with the canonical notion of entropy given above, but 
this result is far from trivial. We see that In Ar, giving equal apriori weight 
to the members of a certain selected sample of quantum states, implements 
the philosophy of the microcanonical ensemble picture. 

In the following sections we are primarily interested in leading-order- 
estimates, that is, estimates of quantities on a scale given by macroscopic 
volumes or by the area of some bounding surface etc. So, by inspecting 



(31) 
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the above formulas, we see that a rough notion which reflects the number of 
"different" states being involved in the reduction of a vector state to a partial 
trace over a certain subvolume, is given by the dimension of the respective 
Hilbert subspace the contributing vectors are lying in. This however needs 
more qualifications. 

In a first step we regard the above representation (formulas (|27p . (|28p of 
the vector $ as a map from the Hilbert spaces H\ , Ti.2 to the corresponding 
subspaces spanned by the (in general non-orthogonal) fa , ipj . These maps 
(more properly, the respective image spaces) are given by C = (cij) and the 
transposed matrix via 

C : e ■ — » c ij e 'j = 0i e t c - (32) 

3 

From linear algebra we know that the dimension of these subspaces is 
given by the rank of the matrices C, C T , with the ranks of C and its trans- 
pose C T being equal. Identifying the matrix C with the corresponding 
abstract map C, it is obvious that this dimension cannot depend on the 
chosen basis. Or put differently, while the subspaces themselves may change 
under a change of basis, the dimension is a unitary invariant. Without tak- 
ing appropriate measures, these subspaces may (in general) have an infinite 
or "too large" dimension. What we are rather implying is the possibility to 
find lower-dimensional essential subspaces. 

Observation/Definition 3.4 The dimension of the subspaces, spanned by 
4>i or ipj as image vectors under the maps C , C T , applied to certain bases 
are equal and independent of the chosen bases in 7i\ , Ti.2- Furthermore the 
logarithm of the dimension is representing a measure of the entanglement 
entropy in leading order under certain favorable conditions. These condi- 
tions may be: All image vectors contribute with roughly the same strength, 
or, rather, their respective weights are not too different. Or, some of the 
matrix elements happen to be very small compared with the rest so that one 
can go over to a purified matrix; in particular, some of the image vectors 
may be very small and can be neglected. In all these cases one will get some 
approximation of the original entanglement- entropy. 

Corollary 3.5 Note that in the Schmidt-basis the dimension of the essen- 
tial subspace can particularly easily be read of. It is the dimension of the 
subspace spanned by the eigenvectors of W\ or W% with non-marginal eigen- 
values Xi . It can of course happen that there is no natural division between 
marginal and non-marginal eigenvalues. This is for example the case when 
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the entanglement entropy is proportional to the volume, a situation we will 
observe in the following. 

In the following sections we will apply this approximative concept of entanglement- 
entropy. 

4 Perturbation Theory of Hamiltonians and En- 
tanglement 

The preceding section shows clearly that the details of the dependence of 
(entanglement-)entropy on volume and/or area as a function of the type 
of global vectorstate are presumably subtle and intricate and need more 
qualifications. This holds the more so if one wants to discuss whole classes 
of models. Therefore, it is a natural idea to study eigenstates of some 
Hamiltonian as described above. In most of the examples we are aware of, 
groundstates of certain physical models have been studied. It hence suggests 
itself to extend this investigation and inspect both ground and excited states 
more closely and try to infer characteristics of the respective entanglement 
entropy from the (general) properties of such eigenstates. 

A general method which suggests itself is perturbation theory of the 
Hamiltonians under discussion. The strategy is the following. We start 
from a Hilbert space, Hy, over a macroscopic volume V and divide it into 
two connected subvolumes, V = V1UV2, both V\ , V2 still being macroscopic 
with 

Hy = Hi® H 2 (33) 

We assume a Hamiltonian, Hy, to be given on Hy of the kind described 
above. This Hamiltonian can be written as 

Hy = H Vl + Hy 2 + H bd (34) 

with Hy i =: Hi the commuting Hamiltonians of the regions V\ , V2 and Hm 
denoting the part of the interaction which comprises lattice points of both 
V\ and V2. That is, we have 

[H 1 ,H 2 }=0, [Hi,H bd }^0 (35) 

It sometimes happens that we have to discriminate between, for example, 
Hi restricted to the subspace Hi and its embedded version, acting on the 
full Hilbert space by tensoring with the unit operator of the volume V2. 
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Definition 4.1 We denote the embedded operators by H\ etc. and the re- 
stricted versions by H\ etc. 

The idea now is to proceed in the following way (we explain it for the 
groundstate). We generally assume that our groundstates are not! degener- 
ate (this is frequently the case). We therefore have the unique groundstate, 
\Po°\ of our reference Hamiltonian, H\ + H 2 , which is the tensor product of 
the groundstates of Hi, i.e. 

*(° ) =4 1) ®4 2) (36) 

One may now hope that one can get information about the structure of the 
groundstate of the full Hamiltonian, H = H\ + H 2 + H bd , by employing per- 
turbation theory with respect to H\ + H2 around the reference groundstate 
*f, 0) . The probl ems and obstacles which show up if one wants to follow this 
route are described below. But anyway, as this strategy may not be entirely 
futile and, to our knowledge, has not been attempted in the past in this 
field, we describe now some of the necessary steps. 

Generically Hi are of "size" Vi (e.g. their respective operator norms or 
(most of their) eigenvalues). That is 

|| Hi \\< Ci-n{Vi) = 0{Vi) , || H bd \\< C'-n(V 1 ,V 2 ) = O (boundary^) (37) 

Therefore one may have the idea to treat Hbd as a relatively small pertur- 
bation of the operators H\ or H 2 . 

The operators H\ , H 2 , defined above over the regions V\ , V 2 , commute. 
In the following we will encounter in various arguments such pieces of the 
total Hamiltonian which commute. Another example is the following. We 
subdivide V\ further into regions V[ , V" and correspondingly for V 2 with 
V\ = V[ U V[. V" is the region in V\ which lies within distance d > 5 {5 
the maximum over the ranges of the various interaction potentials ) of the 
common boundary with V 2 . In certain calculations we choose d macroscopic 
but L S> d S> 5 (where, as usual, we take S as a microscopic quantity). The 
respective Hilbert spaces are H[ , H" with Hi = H[ <8> H". We can now 
define another approximation of the total Hamiltonian H in deleting the 
boundary terms with respect to the interfaces separating V{ , V" on the one 
hand and V 2 ' , V 2 on the other hand, 

H' := H[ + H'{ + H bd + H' 2 + H'> =: H[ + H' bd + H' 2 (38) 

with 

H' bd := H'{ + H bd + H'i (39) 

We now have 
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Observation 4.2 The operators H[ , Hi, , H' 2 commute and 

H = H 1 + fliw(i) + i?6rf(2) (40) 

wii/t the two boundary contributions describing the interaction through the 
interfaces between V[ , V" and V 2 ' , V 2 ■ The advantage is that we now still 
have included the interaction through the common interface between V\ , V2 
in H', the interface we are originally interested in. 

As to such commuting operators we have the following spectral result (cf. 
e.g. [17] or [50]) which goes back to v. Neumann. 

Theorem 4.3 With (for simplicity reasons) H\,H2, ■ ■ ■ commuting bounded 
selfadjoint operators, they all are functions of a common selfadjoint oper- 
ator, A, i.e. Hi = fi(A). It follows in particular that, in case the spectra 
are discrete, it exists a complete set of common eigenvectors for this set of 
commuting Hamiltonians (including multiplicities). 

One problem which however arises immediately if one wants to apply 
perturbation theory of operators is the following (if one is not entirely cav- 
alier as to mathematical rigor). We know from almost every discussion of 
the foundations of statistical mechanics that for macroscopic volumes the 
spectrum of e.g. the corresponding Hamiltonians, while being frequently 
discrete, is nevertheless so extremely dense and/or highly degenerate that 
ordinary perturbation theory is practically useless. A rough estimate yields 
the following qualitative results. The number of eigenvalues (counting degen- 
eracy) of a hermitean matrix is the same as the dimension of the underlying 
Hilbert space. This means in our case 

^(eigenvalues of H v ) = D Ny (41) 

(D the dimension of the local Hilbert spaces TC Xi , Ny the number of sites 
in V). On the other hand, the spectrum of the corresponding Hamiltonian 
extends typically over an interval of order |V| . 

That is, whereas the higher excited states are typically much more de- 
generate and the spectrum is certainly not evenly distributed, a very crude 
estimate yields a typical density of states of the order 0(|y| _1 -D^). This 
prevents the immediate and naive application of ordinary (analytic) pertur- 
bation theory, which works well for perturbations which are small compared 
to the distance of neigboring eigenvalues of the unperturbed Hamiltonian. 
To be more precise, one knows from beautiful results derived by Rellich (see 
[45j ) . and in particular for the finite dimensional case, that for hermitean 
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perturbations the (discrete) eigenvalues and eigenstates are real-holomorphic 
functions in the coupling constant and that this does not only hold for very 
small values (see also [35], [50] or [38]). 

Theorem 4.4 With H £ := Ho + eV selfadjoint for e £ M., Hq , V bounded 
(for simplicity reasons) and H e having purely discrete spectrum, the eigen- 
values Aj(e) and eigenvectors ipi(s) of H £ , with Aj(0), ipi(0) the eigenvalues 
and eigenvectors of Ho, are real analytic functions of e. One can in partic- 
ular choose e = 1. It can however happen that eigenvalues cross (and hence 
degeneracies change in a superficial sense; see the following corollary). 

Corollary 4.5 This has the important consequence that multiplicities be- 
longing to a fixed Aj(e) can only change at e = as analytic functions, being 
identical on a certain interval, are necessarily the same everywhere. The 
other e-values where a singular behavior can occur lie in the complex plane 
away from the real axis. Note however that, as the eigenvalue functions can 
cross at some points, the counting of degeneracies at such points is a matter 
of taste. The dimension of the total eigenspace is of course the dimension of 
the union of the individual eigenspaces belonging to the different Aj(e) which 
meet at that point. 

Remark: The deeper reason why the nasty Puiseux-series can be avoided 
derives from the fact that we have that H £ is self- adjoint for real e (see [45]). 

On the other hand, convergence radii of the corresponding local power 
series expansions happen to be of the order of the distances between the 
points of the spectrum. This prevents to some extent concrete quantitative 
estimates. To see more clearly the true nature of the problem, we can for 
example start from the unperturbed groundstate, \?n , of H^ := Hi + Hi 
and try to infer with the help of perturbation theory the structure of the 
corresponding groundstate of H^ + Hm as a power series expansion 
with respect to the eigenvectors of H(°> . That is, 

*o = E c »-*S ) (42) 

with 

Cm = c$ + c« + . . . , 4 0) = 1 , cgj> = for m ^ (43) 
The first order yields 

c^=V mO /(E^-E^),m^0, and V m0 = (*$\V ■ ^) (44) 
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We see that for perturbation theory to make sense, 

|y m0 |<|^-4 0) l (45) 



Observation 4.6 While is much smaller than Hi or H2 in general, 
it is still a macroscopic perturbation compared to the typically microscopic 
distances between the eigenvalues of Hi . So ordinary perturbation theory is 
not immediately applicable. 

Remark: There exists however a (complicated and tedeous) way to deal with 
such problems to a certain extent (at least in the physics literature) ; see [19] . 

What will however better work is another important method of estimat- 
ing eigenvalues and their behavior under perturbations which does not focus 
so much on the motion of individual eigenvalues under a perturbation but 
rather makes more global and qualitative statements. This method provides 
however no information about the respective eigenvectors, our main point 
of interest. The method is based on the so-called Rayleigh-Ritz-principle 
and/or the Poincare-Courant-Weyl estimates (see [50], [51] or [52]). All 
these statements are based on minimum-maximum- or maximum-minimum- 
estimates and the principle of stronger or weaker constraints on sets of com- 
parison Hilbert space vectors. 

A result, useful in our context, can e.g. be found in [50J, p. 224, called 
the Weyl-Courant-inequalities, which we reformulate here for bounded her- 
mitean operators with discrete and only finitely degenerated spectrum (not 
having zero as an accumulation point). 

Remark: For various reasons the numbering of eigenvalues is different in 
[50J. We start the counting, beginning with the groundstate. 

Theorem 4.7 With A\,Ai operators of the above kind, with sets of eigen- 
values, chosen in increasing order (counting multiplicity), 

El<E\<... ■ El <E\<... (46) 

and 

Eq<E x <... (47) 
the corresponding eigenvalues of A := A\ + A<i, we have the estimates 

E p+q >El+E 2 q , p,q = 0,l,2,... (48) 
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Corollary 4.8 With H -- 
have 

El 

and more generally 

El- || V ||< E p < El+ || V || (50) 
Eq,Eq the groundstates of H,Hq respectively. 

Proof: The lhs of the inequalities follow directly from the theorem and 
\Eg\ <|| V || for all the eigenvalues of V = A 2 . The rhs follows from the 
theorem by interchanging the roles of the operators, that is 

A x = A - A 2 (51) 

and hence 

El +q >E P + E q (-A 2 ) > E p - || A 2 || (52) 

i.e. 

E v < El +q + || A 2 || (53) 

which yields the result by choosing q = 0. □ 

5 The Groundstate of the Hamiltonian 

We begin with the calculation of the entanglement-entropy of the ground 
state, °f the full Hamiltonian over V = V\ U V 2 . In the following we use 
the leading-order identification of entanglement-entropy made in observa- 
tion/definition [33] and the perturbational results of the preceding section. 
In a first step we study the entanglement-entropy of the ground state, ^q, of 
the approximate Hamiltonian, H', introduced in the preceding section (see 
formula (|38p ). We saw that H' can be written as 

H' = H[ + H' 2 +H' hd (54) 

with all the terms on the rhs commuting with each other. Assuming again 
that the ground states are not degenerate we infer from the results of the 
previous section that the ground state energy, E' , of H' can be uniquely 
written as 

E' = E' 0tl + E' 0t2 + E' 0M (55) 

with the rhs the sum of the ground state energies of the terms occurring on 
the rhs of the previous equation. 



Hq + V and V a small perturbation of H$, we 
■ || V ||< Eq < E]+ || V || (49) 
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Remark: Note that the embedded Hamiltonians always have full subspaces, 
belonging to an eigenvalue. For example, H[ has the eigenspace ipQ 1 <g> H'{ 
belonging to the ground state energy E' 1 in the Hilbert space H\ = H[ fS>H'{ 
and ipQ 1 the unique ground state of the restricted acting on H[ ■ 

The need to constantly make these distinctions is a bit nasty and we will be 
a little bit sloppy if no confusion can arise. We then have 

Observation 5.1 The ground state, ^' of H' can now be uniquely repre- 
sented as the tensor product of the ground states of the restricted Hamilto- 
nians, i.e. 

*o = V>o,i ® ^0,2 ® 4>'o,bd (56) 

where 

$, >M €Wi'®W2 (57) 

In order to calculate the partial traces with respect to Hi or H2 we have 
in a first step to develop with respect to a basis of Hi <8> H2 or, what 
amounts to the same, H'i^H'I^H^^H 1 ^- Choosing as bases in the subspaces 
the eigenvectors of the restricted Hamiltonians H[ ^ , H' 2 ^ , H" ^ , H' 2 ^ , 
we can infer the following from the above observation. 

Conclusion 5.2 In the representation of ^/' with respect to the mentioned 
basis in 

n 1 0H 2 = Hi ® n'{ ® H 2 ' H' 2 (58) 

the outermost left and right terms on the rhs remain fixed (no summation). 
The only summation occurs in the boundary term, ip' bd , which is developed 
with respect to a basis in H'[ <8> H 2 ■ Taking for example the eigenvectors of 

H l,(r) ' H 2,(r) We Write 

<Jbd = E <i <i ® €2 (59) 

and 

*o = tfo,i ® (£■••) ® ^0,2 (60) 

We have 

H' bd = H'{ + H'± + H bd (61) 

where the operators occurring on the rhs are all of roughly the same size, i.e. 
of order 0(boundaryi^ 2 ) ■ While H'{ , H 2 commute, the support of over- 
laps both with the support of H" and H 2 and the respective commutators 
are typically different from zero. 
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If we now view as a state over A 1 , the algebra on 7i\ = H[ <8> H" we 

get: 

(%\(A 1)%) = (Vo,i ® < M P ® l)Vo,i ® ^' , bd ) (62) 

where on the lhs 1 is the unit operator on Ti.2, on the rhs it denotes the unit 
operator on Ti!^ . Inserting ip' M = E ■ ^'i <8> ^',2 i n the above expression 
we get 

{%u ® 1)^ ) = J2J2^ c 'M,i ® <ii^ ^o,i ® €1) 

ij Z 

= E M$>,i ® V>o,i ® V£ 1) (63) 

With 6y = Ei^i- 

Conclusion 5.3 The reduced state or density matrix on Hi, corresponding 
to the total vector state, ^' , is 

^ = Ki><<il®< ( 64 ) 
with W" the density matrix on H'{ with matrix elements bij . 

If the local Hilbert spaces have uniform dimension D and with the as- 
sumed finite interaction distance 5, we conclude that the dimension of the 
Hilbert space H'[ is of order 0(D^ bdl ' 2 ^). From the preceding conclusion we 
infer that for the vector ty' Q all perturbations are essentially restricted to the 
boundary region. 

Conclusion 5.4 For lattice Hamiltonians as we have introduced them, the 
groundstate of H' is expected to contain or is scattered over a number of 
eigenstates of Hi of the order 0(D^ bdl - 2 ^). Correspondingly we infer that its 
entanglement- entropy is of order 0{\bd\^\)- 

Remark: The deeper reason why we are able to infer such a general result 
for the approximate Hamiltonian H' is, on the one hand, the sufficient lo- 
calisation of the perturbation in a boundary layer of finite thickness and, 
on the other hand, the uniqueness properties of the groundstate as a ten- 
sor product of the corresponding groundstates of the Hamiltonians of the 
subvolumes. 

Now we come to the groundstate of the full Hamiltonian 

H = H' + + H bd (2) (65) 
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The difference between H and H' is a small perturbation on the scale of H 
or H' as operators, but not on the scale defined by the difference between 
neighboring eigenvalues of H , H' . From the preceding section we know at 
least that 

E' P - A < E p < E' p + A (66) 

with 

A =|| H bd(1) + H bd{2 ) ||= 0(\bondary 1:2 \) (67) 

Our plan is to make an inference from the number of eigenstates of H[ + H" 
which essentially contribute in the representation of ^>' to the correspond- 
ing number which essentially contribute in the representation of ^q, the 
groundstate of the full Hamiltonian. This number should be of the same or- 
der as the corresponding number of eigenstates of Hi, as both sets represent 
complete bases in Hi = H'i <8> H'{ ■ 

We do this in several steps employing the following reasoning. In a first 
step we add the boundary interaction H bd ^ 2 ) to the start Hamiltonian H' 
yielding the intermediate Hamiltonian H[ + H(V" U V 2 ). Its groundstate is 
V"o 1 ® 0o w ^ n 00 groundstate of H(V" U V2). In V\ we have more or 
less the same situation as before with possible perturbations again confined 
(by definition) to the region V". The same argument as before yields an 
entropy for the reduced state over Vi of order 0(\bdi^\)- Now we employ 
the fact that the entropies are necessarily the same on both sides. That is, 
we arrive at 

Observation 5.5 The entropy of the state (j)' reduced to V 2 is of order 
0{\bd 1>2 \). 

Now we employ the localisation properties of -Hm(i) about the interface 
bdi within a small strip of diameter 25, the interface itself having distance 
d^> 5 from the common boundary between Vi and V 2 . From general expe- 
rience, drawn from the foundations of statistical mechanics and many-body- 
theory, we feel allowed to assume that deep inside the region V 2 , i.e. in V^', 
the groundstate, of the full Hamiltonian H should look similar to the 
groundstate, $ of the Hamiltonian H(V" U V 2 ). 

Remark: Note that we always make the assumption (cf. the introduction) 
that our system is not in a quantum- critical state, i.e. correlations do not 
extend to infinity. The latter case would need some extra discussion. 

Concerning the groundstate of the latter Hamiltonian we learned that its 
restriction to V 2 has an entropy of order 0(bd\ 2 )- From the above we again 
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conclude that, in V 2 , ^0 differs from (the groundstate of H') essentially 
in a boundary layer about the interface bd\ 2 - By symmetry we infer the 
same for the region V\ and arrive at 

Conclusion 5.6 The above chain of reasoning leads to the conclusion that 
as a consequence of the spatial localisation properties of H b d , H bd ^ , H bd ( 2 ) 
and certain natural assumptions about clustering or decay of influence and/ or 
interactions, the groundstate of the full Hamiltonian, H , has an entanglement- 
entropy of order 0{bd\ 2 ). 

RemarkiThis general conclusion is corroborated by exact results for (typ- 
ically) low-dimensional models; cf. e.g. the literature mentioned in the 
introduction. 

6 The Higly Excited Eigenstates 

It turns out that for highly excited eigenstates it is more fruitful to adopt a 
completely different strategy, which is strongly inspired by ideas and meth- 
ods taken from statistical mechanics. To simplify the discussion we treat 
the following system. We take a huge box of sidelength L as the total vol- 
ume V. We partition it by a lattice of small boxes, C\, of sidelength I with 
L S> / 3> 5 (5 the range of the interaction in the original Hamiltonian, H. 
I.e., we assume that I is small but still macroscopic; this is the usual assump- 
tion in statistical mechanics. As subvolumes, V\ , V 2 we take certain regions 
in V each of which contains an integer number of such small boxes. I.e., we 
assume (with iVj , iV/ i , iV/ )2 the respective numbers of boxes in V , V.) 

N t = N lA + N h2 , iVj = L 3 /l 3 , |V- 1 = I 3 • iV M (68) 

Each of the small boxes contains Z 3 /a 3 lattice sites of the original lattice. 
We assume of course that the interface, separating V\ and V2, is sufficiently 
regular, i.e. its area is assumed to be of order 0(L 2 ). 

In each of the small boxes, Cj, we take as Hamiltonian, hi, the piece of 
our original total Hamiltonian H with interaction terms confined to Cj, i.e., 
leaving out the interaction terms occurring in H between the different boxes. 
Note that, due to the assumed translation invariance of our interaction, all 
the hi are equivalent as operators. As these small boxes still contain quite 
a few lattice sites, the spectrum of h = hi may still beboth complex and 
degenerated. As new reference Hamiltonians in V , Vj, j = 1, 2, we take 

H' := H' v H'j := ^ K (69) 

V v 3 
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There is now, in contrast to the preceding section, no boundary term 
Hbd, operating in the vicinity of the boundary, bd\2, but the entanglement 
structure may still be quite complex as we will see below. 

Remark: Such truncated systems are frequently discussed and their prop- 
erties exploited in quantum statistical mechanics within the Gibbsian (en- 
semble) approach. See for example [53] , [M] , [55] or [56] . 

As all these hi commute (by construction), the eigenstates and eigenvalues 
of H' or H'- can be built up from the more elementary components belonging 
to the hi. 

So, to begin with, let us start from some macroscopic volume, V , of 
the above kind (V' = V or Vj), the Hamiltonian, H' = Y2v an< ^ some 
eigenvalue, E, sufficiently far away from the groundstate energy Eq = ^2 Ejp 

(i) 

with E- denoting the j-th energy level of the box Hamiltonian hi. We have 

E = Y. E f (70) 

Gi 

(i) 

for certain appropriate combinations of energy levels, E - , of the hi. It is 
here where certain arguments of combinatorial statistical mechanics enter. 

The problem can now be phrased a little bit differently. With ./V boxes 
given and h = hi having the energy levels (counting multiplicity!) E\ < 
Ei < . . . < Ej , . . . , we are interested in the number of ways of distributing 
the energies, Ej, over the N boxes under the constraints 

# = E^ . E = Y. N r E i ( 71 ) 

with Nj the number of boxes having energy Ej. Each such configuration is 
hence characterized by the sequence, (Ni, JV 2 , . . . ,Nj, . . .). We then have 

Observation 6.1 To each fixed configuration (N\,N2, ■ ■ ■ , Nj, . . .) the num- 
ber of ways of distributing the energies Ej over the N boxes under the above 
constraints is 

W=(N\/N 1 \----N j \---) , N = J2Nj, E = Y,NjEj (72) 

From the combinatorics of such expressions one knows that there exists a 
pronounced maximum of W for a special configuration (N\, N2, ■ ■ ■ , Nj, . . .) mt 
(cf. the above cited literature for more details). The constraints can be 
implemented via Lagrange multipliers with, in the end, the multiplier (3, 



25 



belonging to the .©-constraint, turning out to be something like an inverse 
temperature. It is however important (while usually not openly mentioned 
in the literature) that for this and other results to hold, (E — Eq) has to 
be macroscopic, i.e. of order 0(|V|). This implies that, with the individual 
levels, Ej of h, being microscopic, most of the occurring Nj are sufficiently 
large so that Stirlings formula can be applied. After some calculations one 
winds up with the formula (for details see e.g. [M], chapt. 87, p. 364 ff, or 
[53], chapt. 46, p. 150 ff) 

Conclusion 6.2 Under the assumptions being made we have for the most 
probable configuration and the configuration entropy: 

7V J ./iV = e -^/E e ~ /3Ej > ^W max = N - InC^e- 13 ^ + (3 ■ E) (73) 

3 

(with (3 only implicitly given by the first equation). In any case, lnW max 
turns out to be in general proportional to the volume \V\ (N ~ \V\) for 
highly excited states. 

These findings have the following consequences for our entanglement 
problem. With 

*E = J2 c v ( l ,( i El) ®^ E2) > E 1 + E 2 = E (74) 

and (f)!- , wf'^ eigenvectors to the fixed energies E\ , E 2 of H[ , H f 2 , re- 
spectively, this is an eigenvector for H' with energy E. We have just seen 
that the eigenvalues E\ and E 2 are extremely degenerated so that the num- 
ber of terms, occurring in the above sum, can in principle be chosen to be 
of order 0(e v ). 

Note that, in addition, we could also sum over all possible combinations 
of Ei , E 2 with E\-\-E 2 = E but this is not necessary for our argument. We 
can now make various choices. We can for example select a very special and 
simple eigenvector of product type (i.e. all Cjj = except one): 

m E = <ff&) ® ^) , Ei + E 2 = E (75) 

Its entanglement-entropy is of course zero. 

On the other hand, due to the huge degeneracy of all macroscopic energy 
levels of H', we can exploit our above conclusion 16.21 and what we said in 
the preceding sections about our coarse approximation of entanglement- 
entropy. It is easy to choose the Cjj almost evenly spread over the full 
range of degenerated eigenstates belonging to E\ and E 2 in such a way 
that a typical eigenvector to energy E has an entanglement-entropy which 
is proportional to the volume |V|. That is 
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Conclusion 6.3 Due to the huge degeneracy of macroscopically excited en- 
ergy levels of H' , the typical eigenvector, belonging to the class of eigenvec- 
tors of such an energy level, has an entanglement- entropy of order 0(\V\), 
more specifically 

S{W 1 )=a l -\V 1 \=a 2 -\V 2 \=S{W 2 ) (76) 

(a consequence of conclusion WJfy . Furthermore, our preceding discussion 
shows that these states, W\ , W 2 display features we know from statistical 
mechanics. By "typical" we mean, by randomly selecting one of the admis- 
sible eigenvectors from the huge class, we will get such a state with high 
probability as one knows that the combinatorial coefficients, calculated in 
conclusion \6.2[ are extremely large for macroscopic volumes V, Vj . 

Observation 6.4 We remind the reader of our construction of a vector 
state belonging to a canonical equilibrium state of system (1) with the help 
of tensoring with a system (2). Our above findings on higly excited states 
represent, so to speak, the dual version of this observation. Highly excited 
states on V have, as we have seen, a tendency to resemble states on, say, V\ 
which display a marked statistical mechanical behavior (they are of course 
not always true equilibrium states). 

One can now go on and study the full untruncated Hamiltonian, H, 
starting from such a reference Hamiltonian H' . Assuming that by inserting 
the usually very small boundary terms between the blocks Cj and adding the 
boundary Hamiltonian, H^d, nothing spectacular will happen (the ordinary 
assumption in statistical mechanics), the degeneracy of eigenvalues and/or 
the density of states will remain essentially the same. This may be inferred 
from perturbation theory. We therefore arrive at the final conclusion 

Conclusion 6.5 In the case of the full Hamiltonian H , we have essentially 
the same result as in the preceding conclusion for the reference Hamiltonian 
H' . Due to the expected huge density of states and/or the huge degener- 
acy of eigenvalues we have for macroscopically excited eigenvalues, E, that 
a randomly selected (i.e., generic) eigenstate has an entanglement- entropy 
which is proportional to the volume in leading order. 

7 Low-Lying Excited States 

We now discuss the special case that the excited states lie in the vicinity of 
the groundstate, i.e. instead of energy levels fulfilling 

A:=(E-E ) = O(\V\) (77) 
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we deal with excitation energies which are much smaller. It is perhaps 
surprising that in this case we have to use yet another strategy. The same 
general formula 

W = (N\/Ni \ - ■ ■ Nj\ - ■ ■) (78) 

holds of course also in this regime but for example Stirlings approximation 
is no longer applicable as the Nj are in general too small. Even if we 
would ignore this fact (which would presumably only affect the quantitative 
aspects of some estimates) there exists yet another more serious problem. 
The energy constraint (we now denote the energy levels of H' by E^) 

Ej ■ Nj = E' Q + A = N ■ E Q + A (79) 

with A <C \ V\, is more difficult to implement in this regime. It is interesting 
to analyse the consequences of A <C |V|. 

At first glance it seems that we will get the same results as in the previ- 
ous section by applying the same methods (and in the statistical mechanics 
literature known to us we have found almost no remark as to possible prob- 
lems). The method we applied previously is indeed very general but there 
exists a subtle point. The Lagrange multiplier is only implicitly defined 
via the constraint 

E'/N = J2 E i- e ~ PE > I £ e ' PEj ( 8 °) 

i.e., it regulates the average energy per box Hamiltonian, hi, in form of a 
canonical distribution over the energy levels of h. The Ej are in general 
not known in detail but one may infer that with (E' — E' ) = 0(\V\) both 
sides are of the same order for finite (3 so that it is reasonable that we 
can find some definite value for which the implicit equation for (3 can be 
fulfilled. But we now have (£" — E' )/N <C 1 and we conjecture that in this 
regime the above implicit equation can only be fulfilled for (3 3> 1 or (3 — ► oo 
(which seems to be quite natural, given the obvious similarities to statistical 
mechanics. A "thermal" state near the groundstate has by definition a low 
temperature). 

Observation 7.1 For (E' — E' Q ) small, i.e. (E' — E' ) <C 0(|V|), the param- 
eter (3 becomes very large. For these values it becomes difficult to reliably 
estimate the terms in the occurring variational equations which are now com- 
binations of very large and very small terms. Note in particular that for (3 
large 

Nj = N-e-P E i/J2 e ~ PEj ( 81 ) 
becomes very small compared to N. 
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Therefore we choose another strategy which is better adapted to this 
situation. We simply catalogue the low-lying excitations of H' directly, 
beginning with the groundstate. We have 

Observation 7.2 1) For the groundstate we have: 

E' = N ■ Eq , no degeneration (82) 

2) For one box hamiltonian excited we have: 

E\ = (N- l)Eo + Ei , W{E[) = N (83) 

3) Two levels excited; there are two possibilities, E-i = Ej or Ei ^ Ej. We 
have 

E' u = (N — 2)E + 2E t , W(E' ri ) = (N ■ (N - l)/2) (84) 

or 

E\j = E t + Ej + (N — 2)E , W^) = N(N - 1) (85) 

etc. 

Remark: Note that these results of course coincide with the general formula, 
if we insert the corresponding iVj. 

We see the following. Already for the lowest excited levels of H' we 
have a degeneracy, W(E' i ) = 0(\V\) or a small power of \V\. Repeating our 
previous arguments we infer 

Conclusion 7.3 Already the lowest excited levels of H' have a degeneracy 
of order 0(\V\) or a small power of\V\, entailing that we can construct cor- 
responding eigenstates having an entanglement- entropy of order 0(ln|V|). 
Note that H' does not contain the boundary term H^. For the full H we 
hence expect a generic entanglement- entropy of order 0(ln \ V\ ■ |6di2|) where 
the additional factor \ bd\2 \ comes from the term as this term leads to a 
further splitting proportional to the area of the boundary as was for example 
the case for the groundstate. 

The situation changes slightly if we go over to higher excited levels. For, 
say, k levels excited the two extreme cases are: 1) all k levels identical or, 
2) all levels being different. The intermediate class comprises cases where 
some of the Ej coincide. We have the following estimate 
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Corollary 7.4 // k levels are excited, , . . . , Ei k , with repetitions allowed, 
we have the following estimate with respect to the degeneracies of the low- 
lying eigenvalues of H' 

(N\/k\ ■ (N - k)\) < W k < {N ■ (N - 1) • • • (N - k + 1)) (86) 

For N still very large compared to k, this entails 

W k = 0{N k ) and In W k = 0(k • In N) = 0(k • In \V\) (87) 

For the full Hamiltonian H we again get an additional multiplicative factor 
of order 0(| 6cZi2 1 ) - 

8 Appendix 

A Infinite-Dimensional Local Hilbert-spaces 

We assume that the spectra of the local and global hamiltonians are discrete 
but are now not necessarily bounded from above. We still assume that Hy 
is bounded from below (existence of a groundstate!). So let H be such a 
Hamiltonian on a separable Hilbert space, TL. We select a certain (countable) 
basis, ej, and choose certain subspaces, 7i n , spanned by the basis vectors, 
ei, . . . , e n . The projector on 7i n is denoted by P n . Then 

P n HP n (88) 

is a bounded operator on 7i n . Being a little bit more general, if Hy is the 
Hamiltonian on Hy := 7ii 1 (8> ■ ■ ■ <8> 7~li N , we select in each Tti v the subspace 
Ti\ n ^ , spanned by e±\. . . , e$ . From these local pieces we compose the 
subspace, H n , in Hy, i.e. 

H n := W£° <g> • • • <g> Hj} (89) 

and denote again the projector on this subspace by P n . In the same manner 
we take 

ff (n) :=P n -HyP n (90) 

as finite dimensional Hamiltonian on these 7i n C Hy. 

Now we have to discuss what happens if we take the limits 

H n -► Hy , ff (n) -»• Hy etc. (91) 
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Note that by construction all the are now bounded, even finite dimen- 
sional, but the limit Hamiltonians are in general unbounded. We do not 
want to be too tedious concerning technical details of functional analysis at 
this point. Suffice it to say that a reasonable concept of operator conver- 
gence in this context is convergence in the resolvent sense, i.e., instead of 
dealing with unbounded operators we deal with their bounded resolvents 

(H - z)- 1 , Im(z) ^ (92) 

Under quite weak assumptions, which are in general fulfilled in our con- 
text, (strong) resolvent convergence can be assumed (cf. [37], sect. VIII. 7). 
As a consequence we have the following result (discrete spectrum): 

Observation A.l With H having discrete spectrum and if H n — > H in 

strong resolvent sense, we can find to each eigenvalue E of H an interval 
(a,b) so that E is the only spectral value of H in (a,b). With Pf^b) ^ e 
spectral projections of H n on the interval (a, b) we have 

P (al)^ ~+ P M*P = Pe^P (93) 

for all if) G H. We get even stronger results if we assume convergence in 
norm-resolvent sense. 

We now sketch how results about entanglement-entropy, derived for large 
but nevertheless finite dimensional systems, could be transferred to the gen- 
eral case. But as there are several quite delicate technical steps involved, 
which to rigorously prove would need quite an amount of mathematical in- 
put, we refrain from giving all the intricate mathematical details at the 
moment. From a physical point of view the strategy seems to be quite 
reasonable. 

In a first step we have to guarantee that for example in the case of 
groundstates 

- * (94) 

In the following sections we regard these vectorstates as states on a restricted 
region, V\, and the corresponding Hilbert space or observable algebra. We 
then have for the respective density matrices over 7i\ 

w[ n) -> Wi (95) 

in the form 

Tr{w[ n) ■ A) -> Tr(Wi ■ A) (96) 
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What we need is a result like 

W[ n) ■ In w[ n) -»• Wi ■ In Wi (97) 
in a suitable topology so that we may get in the end 

Tr(w[ n) ■ In W[ n) ) -> 2V(Wi • In Wi) (98) 

that is 

S(w[ n) ) -> S'(Wi) (99) 

Remark: Some of the necessary technical arsenal can be found for example 
in [38], see also the seminal paper by Wehrl (|39j), in particular section IID 
about continuity properties of entropy as a function defined over the density 
matrices, or the more recent |40| . 

Note that in general the entropy is not continuous if the density matrices 
are equipped with the trace-norm 

|| Pi ~ pi \\tr-= tr(\pi - p 2 |) (100) 

It is an important observation that this property can be restored if one 
works with the subset of density matrices having finite energy. This means, 
one chooses a fixed reasonable Hamiltonian, H, so that the canonical Gibbs 
state, ap, at some inverse temperature, /3, has finite trace and finite mean 
energy 

tr{a p -H) = E (101) 
The subset of admissible density matrices is then given by the condition 

tr(p-H)<E (102) 

This is a reasonable condition for e.g. entanglement entropy. Taking for the 
bipartite system a Hamiltonian of the form H = H1 + H2 (i.e. no interaction 
across the common boundary) with reasonable Hi, it suffices that the Hi are 
in the domain of definition of the pure state ^ as we then have 

tr(Wi ■ Hi) = (*, (Hi ® 1)10 = E 1 (103) 

for some E' . So, adapting (3, these particular density matrices always lie in 
appropriate subsets. 
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B Proof of Theorem ED] 



Viewing the matrices C := (cj m ) , C* := (C as operators from TC2 — ► 
Hi , Hi — > H2 respectively, we have 

Wx = C ■ C* , W 2 = C* -C (104) 

If ipi is an eigenvector of W2 = C* ■ C, i.e. 

C* ■ C o ^ = A, • Vi (105) 

it follows that 

CC*C o ^ = Ai • C o V (106) 

i.e., Aj is eigenvalue of W% = CC* with eigenvector Ctpi and vice versa. 
Furthermore, if Aj 7^ 0, the degeneracy is the same with respect to C*C and 
CC* . This follows from 

(C#|CV?) = (^IC*^ 2 ) = Ai(^IVl) (107) 

that is, with t/^ 1 , ipf orthogonal eigenvectors to the eigenvalue Aj / 0, 
Cipj , Cw are also orthogonal and non-vanishing and the same holds in 
the other direction. 

By the same token we infer from the positive definiteness of C*C and 
CC* that all eigenvalues are > 0. The normalisation of the vector $ as a 
state on A and A4 implies that the trace norm of Wi is one. This proves the 
theorem. 

Acknowledgement: I thank the referees for various helpful comments. 
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